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^ Abstract 

Supersymmetry plays a fundamental role in the radiative stability of many inflationary models. 

^ Spontaneous breaking of the symmetry inevitably leads to fields with masses of order the Hubble 

scale during inflation. When these fields couple to the inflaton they produce a unique signature 

in the squeezed limit of the three-point function of primordial curvature perturbations. In this 

paper, we make this connection between naturalness, supersymmetry, Hubble-mass degrees of 

freedom and the squeezed limit precise. To study the physics in a model-insensitive way, we 

develop a supersymmetric effective theory of inflation. We use the effective theory to classify all 

possible interactions between the inflaton and the additional fields, and determine which ones 

naturally allow large non-Gaussianities when protected by supersymmetry. Finally, we discuss 

J^j the tantalizing prospect of using cosmological observations as a probe of supersymmetry. 
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1 Introduction 



Supersymmetry (SUSY) [1] is a powerful way to protect scalar fields from quantum corrections. If 
the symmetry is unbroken, it controls dangerous radiative effects by enforcing exact cancellations 
between boson and fermion loops. Even if the supersymmetry is broken at low energies — as it 
has to be if it is to describe the real world — the appearance of supersymmetry at high energies 
can still help to regulate loop effects. In inflation [2], SUSY is often required to achieve techni- 
cal naturalness [3] of the quantum-corrected inflaton action [4, 5]. Hence, even if SUSY is not 
discovered at the weak scale, naturalness motivates SUSY in inflation. Can we use cosmologi- 
cal experiments to probe supersymmetry even if it is out of reach of particle colliders? What 
are generic signatures of supersymmetry during inflation? In this paper we will address these 
questions. 

An inevitable consequence of SUSY during inflation is the presence of additional scalar fields, 
with masses of order the Hubble scale H. The appearance of the Hubble scale as a preferred 
mass scale is robust and simply related to fact that the size of spontaneous SUSY breaking during 
inflation is determined by the vacuum energy, (J-x) 2 = 3M 2 \H 2 . Gravity mediates this breaking 
to the scalar sector, leading to the relation 

m ~ &L rsj H . (1.1) 

M p i 1 > 

This effect is familiar from the supergravity (SUGRA) eta problem [6], where the induced mass 
for the inflaton field, m? ~ H 2 , threatens to end inflation prematurely. Inflation will only last 
for at least 60 e- folds if the inflaton mass is tuned to be smaller than H, or if a smaller mass is 
protected by an additional global symmetry. This makes it clear that additional fields without 
protective global symmetries will inherit Hubble scale masses from SUSY breaking. It is less 
clear that these fields can have observational consequences. 

We are looking for an observable that is sensitive to the number of light degrees of freedom 
during inflation, and that, ideally, would probe the masses of the additional fields. A promising 
candidate is the squeezed limit of the three-point function of primordial curvature perturbations £ : 




k s 



This describes the correlation between a long-wavelength fluctuation C,l = C(&l) an d two short- 
wavelength fluctuations Cs = C(^s)> i- e - Cl modifies the short-scale power spectrum (Cf)- How- 
ever, in single-field inflation, constant Cl is a pure gauge mode and hence isn't locally observable. 
The long-wavelength mode therefore has no effect on the short-wavelength power in the zero 
momentum limit ki — > 0. Any physical effects are proportional to {Cl, V 2 Cl} oc k\C,L and hence 
suppressed by (^l/^s) 2 [7-9]. On the other hand, if the fluctuations we observe were produced 
by a second light field, the superhorizon fluctuations of the second field are locally observable 
at the time when they are converted to curvature fluctuations. Large non-Gaussianity is then 
possible in the squeezed limit. These facts are summarized by the momentum dependence of the 
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three-point function in the squeezed limit 

fc lim(C kL C kSl Ck S2 )oc^. (1.2) 

Here, single-field inflation corresponds to a = 1 [9], while multiple light fields with m <C H can 
lead to non-Gaussianity with a = 3 [10]. 1 Until recently, no example was known with a different 
from 1 or 3. We see that measurements of the squeezed limit have the potential to determine 
whether one or more fundamental fields were relevant for producing the primordial fluctuations. 

How does the presence of Hubble-mass degrees of freedom during inflation affect the squeezed 
limit? For the massive 'isocurvaton' fields a to have an impact on observables, their fluctuations 
have to be converted into curvature fluctuations £ — either during inflation or afterwards. Re- 
cently, Chen and Wang [11] observed that direct couplings between the inflaton <p and a massive 
isocurvaton a lead to a squeezed limit with the following momentum dependence: 

3 [9 m% 

a= 2 + V4-#- (L3) 

These models of quasi- single- field inflation (QSFI) therefore allow for a scaling in the squeezed 
limit that is intermediate 2 between the scaling of single-field and multi-field inflation: | < a < 3. 
The appearance of the mass m a in eq. (1.3) has an intuitive explanation: massive fields decay 
when their wavelengths becomes larger than the Hubble radius. At the time when the short 
modes cross the horizon, the amplitude of the long mode is therefore suppressed relative to the 
amplitude when it crossed the horizon. This affects the momentum dependence of the squeezed 
limit by an amount that depends on m a (see §4.1 for a derivation). 

The squeezed limit is both theoretical clean and observationally relevant. Future cosmic mi- 
crowave background (CMB) data [15] will span a large range of scales and therefore provide access 
to the squeezed limit of the primordial perturbations. These measurements will be complemented 
by large-scale structure (LSS) observations [16]. In particular, in recent years, scale- dependent 
halo bias 3 [17, 18] has emerged as a sensitive probe of primordial non-Gaussianity [19, 20]. In 
this case, the signal is dominated by the squeezed limit of the three-point function [21]: 

Ab(k) oc ^ . (1.4) 

The combination of CMB and LSS data therefore suggests the exciting possibility of using obser- 
vations to probe Hubble-mass degrees of freedom during inflation. As we explained, these fields 
have their natural home in supersymmetry. In this paper, we will make this relation concrete. 
Our results will solidify the connection between naturalness, supersymmetry and the squeezed 

1 Here, and in the rest of the paper, we are assuming perfectly scale-invariant fluctuations, fc 3 (0?) oc k° . It is 
straightforward to restore (percent-level) deviations from scale-invariance, fe 3 (^|) oc fc ns_1 , in all of our results. 

2 Such intermediate scalings may also be engineered in some curvaton scenarios [12, 13], as well as in theories 
with non-Bunch-Davies initial states [14]. 

s Galaxies reside in dark matter halos. Halos are biased tracers of the underlying dark matter density field: 
5h = b<5 m . The halo bias b(k) = (<5 m <5h}/(<5m<5m) develops a characteristic scale-dependence in the presence of non- 
Gaussianity with a non-trivial squeezed limit of the primordial three-point function. We can hope to detect this 
effect by measuring the scale-dependence of the galaxy two-point correlation function on large scales. 
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limit of non-Gaussianities. To use cosmological observations to find evidence for a new spacetime 
symmetry of Nature is, of course, a very tantalizing prospect. 

We have two primary goals in this work: 

1) We wish to develop a broadly applicable framework for incorporating supersymmetry in 
effective theories of inflation. The starting point for this investigation will be the effective 
theory of inflationary fluctuations of Cheung et al. [22] — a theory of Goldstone bosons 
of spontaneously broken time translations. We will show how the Goldstone fields are 
embedded in a supersymmetric theory (see also [23, 24]). Along the way we will have to 
understand how to couple higher-derivative theories to supergravity [25], how to cancel 
tadpoles, and how to characterize the couplings of matter fields to auxiliary supergravity 
fields. The end product will be the supersymmetric effective theory of inflation. This is an 
effective theory of the inflationary fluctuations and therefore complementary to the large 
body of work on SUSY theories for inflationary backgrounds (see e.g. [5, 26]). 

2) We then use this theory for a model-independent description of the interactions between 
the inflaton and additional Hubble- mass fields. The supersymmetric framework allows us 
to analyze naturalness in a UV-insensitive way. One of our main goals will be to classify 
all technically natural SUSY implementations of quasi-single-field inflation. 

The outline of the paper is as follows: In Section 2, we construct a supersymmetric general- 
ization of the effective theory of inflation of Cheung et al. [22]. We will illustrate the formalism 
with two important examples: in Section 3, we derive supersymmetric models of inflation with 
small sound speed, while in Section 4, we present the first explicit models of QSFI in SUSY. 
We classify all possible variations of QSFI and determine which scenarios can lead to naturally 
large non-Gaussianity when protected by SUSY. We discuss the observational signatures of these 
theories. Sections 3 and 4 are self-contained and can therefore be read separately and in any 
order. We conclude with Section 5. 

Two appendices supplement the computations of the main text: In Appendix A, we show 
how to couple higher-derivative theories consistently to supergravity. This is a largely non- 
technical summary of the results of a companion paper [25]. We use results from this appendix 
in Section 3. In Appendix B, we develop simple estimation techniques to derive the amplitude 
and the squeezed limit of the bispectrum for models of quasi-single-field inflation. These results 
are required to reproduce the results of Section 4. 
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2 The Supersymmetric Effective Theory of Inflation 



2.1 Adiabatic Fluctuations as Goldstone Bosons 

We begin with a quick review of the effective theory of inflation of Cheung et al. [22] (see also [27]). 
Readers familiar with this work may skip directly to the next subsection where we generalize the 
theory to include supersymmetry. 

The first step in constructing effective field theories (EFT) is identifying the relevant degrees 
of freedom for the measurements of interest. Here, we are interested in an EFT that describes 
CMB observables. We assume that the origin of the CMB temperatures fluctuations can be 
traced back to quantum fluctuations during inflation that froze and became classical when their 



frequencies uj matched the Hubble expansion rate H. 
inflationary fluctuations valid at energies oj ~ H. 

To construct this EFT, we start from the crucial 
insight that time-dependent FRW backgrounds, such 
as the inflationary quasi-de Sitter spacetimes, sponta- 
neously break time-translation invariance. The relevant 
degree of freedom of the EFT is then the Goldstone 
boson associated with the symmetry breaking. We in- 
troduce the Goldstone mode as a spacetime-dependent 
transformation along the direction of the broken sym- 
metry, i.e. we perform a spacetime-dependent time shift 



The goal is therefore an EFT for the 




cp = t + 7r(x) , 



(2.1) 



where t — > t + £o(x) and tt — > tt — £,o( x )- We note that the field tt parameterizes adiabatic 
perturbations, i.e. perturbations corresponding to a common, local shift in time for a set of 
matter fields 4> a , a = 1, • • • , N, 



f>a(x) = 0o(* + T(x)) - 0«(t) • 



(2.2) 



Moreover, the field tt is related to the comoving curvature perturbation £ by a gauge transfor- 
mation, 

C = — Hit . (2.3) 
The low-energy expansion of the field (p(x) = t + n(x) is the effective theory of inflation [22] 



(2.4) 



Slow-roll inflation. At lowest order in derivatives, we have a (time-dependent) potential 
and a kinetic term for (p, 

C cS = -A(<p) + c( V )(3^) 2 + • • • (2.5) 

To cancel tadpoles (i.e. terms linear in the fluctuations tt), the coefficients in eq. (2.5) are fixed 
in terms of the Hubble parameter: A(ip) = M^(3H 2 (tp) + H) and c((p) = M^H. This ensures 
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that we are expanding around the correct inflationary background H{t). Eq. (2.5) is, of course, 
nothing but slow-roll inflation in disguise 

£ s ,. = M 2 l H(d II <p) 2 -M 2 l (3H 2 ( V ) + H) <- -\{d^f - V{4>) . (2.6) 

A drastic simplification occurs in the effective action when we take the so-called decoupling limit. 
In this limit we ignore the mixing between the matter fluctuations n and the metric perturbations 
5g^ v — i.e. we evaluate the Goldstone action in the homogeneous background spacetime: g^ v — > 
g^ u . This approximation only leads to slow-roll suppressed errors in the correlation functions 
evaluated at horizon-crossing u ~ H. Since the effective theory of inflation is most powerful 
when it describes large non-Gaussianities arising from matter self- interactions, we are justified 
to work in the decoupling limit — e.g. 

(d^) 2 = g^d^dvtp -> -1 - 2tt + (d^ f . (2.7) 

The action (2.6) then becomes 

£ s . r . = M^H^irf . (2.8) 
We find that it is massless and purely Gaussian for slow-roll inflation in the decoupling limit. 

Theories with small sound speed. To go beyond the simple free-field theory, we add higher- 
derivative terms. For single derivatives acting on ip = t + it, the leading deformation of the 
slow-roll action is 

C Cs = £ 8 . r . + \M%{ip) [{d^f + 1] 2 . (2.9) 

Notice that we added '+1' to the operator in eq. (2.7) before squaring. This was necessary in 
order to avoid reintroducing a tadpole for it. In other words, without loss of generality, we chose 
the cancel the tadpole once and for all at lowest order and only add new operators without 
tadpoles. Written in terms of the Goldstone n, we see that the operator in eq. (2.9) modifies the 
kinetic term, but not the gradient term 

C Cs = —{M 2 X H - 2M 2 V + M^HidiTr) 2 + ■■■ (2.10) 

This induces a sound speed for the propagation of n, 

L = l- 2M l (211) 

with a large value for the parameter M2 corresponding to a small sound speed c s < 1. Substitut- 
ing (2.7) into (2.9), we see immediately that a small sound speed relates to large interactions [22] 

C Cs = --- + 2Ml(ir 2 -\ir{d^) 2 + ■■■) . (2.12) 

This non-linearly realized symmetry is particularly clear in the EFT approach. 

P(X) theories. The generalization of the single-derivative Goldstone action to so-called 
P{X) theories [28-30] is straightforward: 

00 

£p(x) = E h M n&) t(^) 2 + If • (2-13) 

n=0 
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Here, the operators proportional to M n start at order n in ir and end at order 2n. Only two in- 
dependent operators — proportional to M2 and M3 — therefore contribute to the cubic Lagrangian 
for 7T. 

Higher- derivative theories, ghosts and galileons. The number of operators increases rapidly 
if we allow for the possibility of more than one derivative acting on 99. Higher derivatives lead to 
interesting theories such as ghost inflation [31] and galileon inflation [32]. We won't treat them 
in this paper, although our formalism of course applies to all of these cases. 

Multi-field generalizations. Finally, we could consider coupling the adiabatic mode ip to ad- 
ditional (isocurvature) fields a. To keep these extra fields naturally light either requires additional 
global symmetries or supersymmetry (with a modest amount of fine-tuning) [23]. 

Super symmetry. In the rest of this section we will extend the effective theory of inflation 
to include supersymmetry. In §2.2, we introduce the basic ingredients of the theory and explain 
our assumptions. In §2.3, we present a systematic treatment of effective actions with global 
supersymmetry. We discuss supergravity couplings in §2.4. Finally, in §2.5, we comment on 
models with additional matter fields. 

2.2 Supersymmetry Breaking during Inflation 

It will be crucial for the construction of the SUSY EFT to understand how SUSY is broken 
during inflation. 

Vacuum energy. The primary source of spontaneous SUSY breaking during inflation is, 
of course, the positive vacuum energy associated with the de Sitter background. Although the 
vacuum energy is present in any model of inflation, it is only by introducing supersymmetry that 
it breaks a symmetry. We characterize this effect through a SUSY breaking spurion 

x = x + V2e^ x + e 2 F x , (2.14) 

with nearly constant F-term 4 , |.Fx-| 2 ~ 3M 2 ^H 2 . The action for X in the rigid limit is 

C x = j d 4 9X j X- (^j d 2 9V3M p iH X + h.c.^ . (2.15) 

As written, this action does not generate a mass for x. One could add higher-order terms in 
the Kahler potential — e.g. 5K = A~ 2 (X^X) 2 — to stabilize the vev at x = 0. In any model 
with m x ^> H, one can arrive at the same results without changing the action by imposing the 
constraint X 2 = [35]. In the following, we will assume that the classical vev of x vanishes. 

For a given inflationary model, the inflaton may or may not be part of the multiplet X. 
However, supersymmetry is most useful when it protects the mass of the inflaton to the lowest 
possible scale, namely the Hubble scale H. The F-term for X breaks SUSY at a much higher 
scale, making it more difficult to protect the inflaton mass when the inflaton is in the same 

4 We do not consider D-term breaking. In fact, mostly D-term breaking has been shown to be incompatible with 
well-defined current multiplets [33, 34]. By restricting to .F-term breaking, we are therefore only ignoring models 
with comparable D-term and F-term contributions. 
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multiplet (see [36, 37] for future discussion and examples). We will therefore make the further 
assumption that the inflaton <j> is in a separate multiplet, 

$ = <t> + V26ij; 4> + e 2 F , (2.16) 

where 

<t> = -^(a + i<p c ) , with ip 2 c = 2M* i \H\{t + ir) 2 . (2.17) 

We note that supersymmetry inevitably adds a new real scalar degree of freedom a. The inter- 
actions between a and it will be of considerable interest for the phenomenology of inflationary 
fluctuations. 

Time- dependent backgrounds. The feature that distinguishes inflation from pure de Sitter is 
that time translations are spontaneously broken by the inflationary background. This necessarily 
leads to a second source of SUSY breaking. We can see this directly from the SUSY algebra 
{QaiQa} = ^aa^ 3 ^- ^he breaking of time translations means that the generator Pq is broken. 
To be consistent, supersymmetry must also be broken. This can also be understood from the usual 
argument that a supersymmetric state must have zero energy. A rolling scalar field generates 
positive kinetic energy and breaks SUSY even in the absence of potential energy. Finally, we note 
that this type of SUSY breaking is reflected in the supersymmetry transformation of the spinor 
field 

= iVWtdrf + V2^F . (2.18) 

The fermion transforms inhomogeneously and, hence, supersymmetry is spontaneously broken 
by the time-dependent inflaton vev, ((f)) ^ 0. In superspace, the breaking manifests itself as 
a non-zero 00 component of $, i.e. ($)\ gS = ia^d^) = -a°(M^\H\) 1 ^ . This type of SUSY 
breaking is not usually considered in particle physics since it breaks Lorentz invariance. We will 
refer to this type of SUSY breaking as 'Lorentz breaking'. 

2.3 Supersymmetric Actions 

In a rigid theory in flat space, we can always decouple the actions for X and <J>, 5 

£ = . (2.19) 

What are the constraints on the form of the Lagrangian £$? Having organized the multiplet 
in terms of tp = t + ir, we have ensured that $ is invariant under local time translations. As a 
result, any supersymmetric action written as a function of <3? and & is time translation invariant, 
including terms with any number of derivatives. However, as we saw in §2.1, tadpole cancellation 
imposes additional constraints on the form of the action. This is where the EFT of inflation differs 
from the effective action for Goldstone bosons of global symmetries, where the only constraint 
on the action is that it is invariant under the symmetry. Finally, we will impose an additional 

Of course, one cannot decouple X and t completely, otherwise inflation would never end. However, this 
coupling may arise indirectly through couplings to fields that are sufficiently massive during inflation that they 
can be integrated out, leaving the decoupled action to good approximation. These fields may become light or 
tachyonic in order to end of inflation without affecting the fluctuations during the e-folds that are probed by the 
CMB and LSS. 
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shift symmetry, ir — > ir + const., on the action to ensure at least 60 e-folds of scale- invariant 
fluctuations for £ = —Hit. The most general 6 Lagrangian of $ with this symmetry is 

U = J d 4 eK($ + &,d^,d^,---) . (2.20) 

For example, in slow-roll inflation in the decoupling limit, the unique action is 

£ s .r. = J d i 9\{$> + $t)2 + Cx . (2.21) 

More general theories are characterized by higher-derivative corrections to (2.21). For instance, 
in §3, we will study the following example 

£c s = [d A 0l($ + &) 2 \c x + — d^d»&] + C x , (2.22) 
J M P i\ H \ J 

corresponding to a supersymmetric theory with non-trivial sound speed. 
2.4 Coupling to Supergravity 

To describe the dynamics of inflation, we must couple the theory to gravity. For a supersymmetric 
theory, this means coupling to supergravity. The decoupling between <3? and X that was possible 
for the rigid theory in flat space is no longer possible. 8 

Tadpole cancellation. The decoupling between 3> and X is broken by the time evolution 
of the vacuum energy, H ^ 0. In order to reproduce the correct time-dependent vacuum energy, 
A(i£>) = Mpj(3ff 2 ((/?) + H), we include a small superpotential coupling between <3? and X, 

C x ^ Jd 4 6X j X - (^j d 2 ev{-i<5>)X + h.c^j , (2.23) 

where |f(c/?)| 2 = M^(3H 2 (ip) + H). This is equivalent to including couplings that cancel the 
tadpoles for ir at all values of t. Additional couplings may be required to cancel tadpoles that 
may arise for the inflaton partner a. 

Curvature couplings. Further couplings between $ and X generically arise in supergravity 
3 aKo 2 + --- + f3 [ d 4 9-^($ + &) 2 X^X + ■■■ , (2.24) 



6 We have actually imposed a slightly stronger condition, namely that this is a symmetry of the Kahler potential 
itself. Here, we are anticipating that supergravity effects can break symmetries of the action that are not symmetries 
of the Kahler potential. 

7 fn flat space, there is no need for additional terms to cancel tadpoles for ir. Because of the shift symmetry, 
any terms linear in n are total derivatives in fiat space. However, we will be most interested in FRW backgrounds, 
where tadpole cancelation requires that we consistently couple to gravity. 

8 It is amusing to note the similarity between the SUSY EFT of inflation and Goldstini [38], as both describe 
two independent SUSY breaking sectors that are coupled only through gravity. These theories are different in two 
respects: First of all, in inflation, the massive scalar a will not, in general, be sufficiently massive that it can be 
integrated out and will play a critical role in the phenomenology of the EFT. Second of all, one of the sources of 
breaking is 'Lorentz breaking' rather than F-term breaking. Even in models in which a is sufficient heavy that it 
can be integrated out, the Goldstino multiplet takes the form <1?nl = ipo-°ijj/Al + tp + ... , where i/j is the Goldstino 
of the Lorentz breaking sector and Aj, is the scale where the symmetry is broken. Unfortunately, since fermions 
don't freeze out, the Goldstini are unobservable. 
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where a and (3 are dimensionless coupling constants of order unity. Notice that only the non- 
shift-symmetric part of <5, i.e. a oc <3? + <J?t, receives corrections, while the potential for the 
shift-symmetric part, i.e. ip oc remains protected. Although curvature couplings, such as 

the first term in (2.24), do not couple the fields directly, Einstein's equations relate the spacetime 
curvature 1Z to the vacuum energy generated by X. This indirectly couples X and $ + The 
second term in (2.24) represents Planck-suppressed couplings between the two sectors. Using 
Tx ~ V3M p \H and 1Z = —12H 2 (here we are using the sign conventions of [39]), we see that 
both terms in (2.24) contribute a mass to a of order H. 9 

Auxiliary fields in supergravity. Minimal supergravity has two auxiliary fields, a complex 
scalar M and a real vector b^ [39] (see Appendix A). These fields can also couple to the inflaton. 
When M and acquire vevs, they can, in principle, affect the inflationary dynamics. For models 
of particle physics (e.g. gravity or anomaly mediation), the contribution to the soft masses from 
M is extremely important: in this case, a vev of the superpotential W leads to a universal 
contribution of the form (M) ~ (W)/M^. Such a vev, in fact, has to be turned on in order 
to achieve a small value for the present day cosmological constant, i.e. Wq ~ TxM p \ and hence 
(M) ~ Tx/M p \. Therefore, SUSY breaking effects communicated by M are comparable to those 
from Planck-suppressed mixing between sectors. In inflation, the situation is very different: The 
vacuum energy plays a crucial role and should not be cancelled by the vev of the superpotential. 
As a result, one requires that (W) <C TxM p \ and therefore (M) <C Tx /M p \ ~ H. This makes the 
contributions from M subdominant relative to the generic curvature couplings from supergravity. 

In fact, under the assumptions stated in §2.2, we find that the vevs of the auxiliary fields M 
and F are strictly zero. The basic reason for this is as follows: The action for the SUSY effective 
theory of inflation has an i?-symmetry. Under this symmetry, X has -R-charge 2 and $ has R- 
charge 0. As a result, both <p and Tx have charge zero and neither vev breaks the .R-symmetry. 
However, both M and F are charged and a vev for either would break the symmetry. Since M 
and F are auxiliary fields, any vev should be proportional to an i?-symmetry breaking parameter. 
The only available quantity is the vev of x. However, if (x) = 0, as we assume throughout, then 
(M) = (F) = 0. 

In contrast, there is no general principle that requires (b^) = 0. However, in all of the models 
that we will discuss in this paper, we will actually find this to be the case. More generally, a 
vev for b^ breaks diffeomorphism invariance. Because time diffeomorphisms are broken, bo can 
acquire a vev. A vev for bo also breaks SUSY and must be related to the "Lorentz -breaking" 
of supersymmetry. Because this is a subdominant contribution to the vacuum energy, (&o) "C H 
during inflation. 

Supergravity for effective theories. In this section we have merely outlined the contributions 
to the action that may arise from supergravity. Determining the precise form of these corrections 
requires a complete treatment of supergravity, including a careful treatment of higher-derivative 

9 In the limit, M p \ — > oo, the second class of corrections naively seem to vanish, while the first class remains. 
However, for the background to solve Einstein's equations for any finite value of M p i, we must take Tx — > oo in 
the same limit. For this reason, both types of terms survive the decoupling limit and contribute terms of order 
H 2 . This non-decoupling of SUSY breaking is necessary. In de Sitter space, the transformation of the gravitino is 
required for any action to be supersymmetric and therefore de Sitter space must break supersymmetry explicitly 
in the decoupling limit. 
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contributions to the Kahler potential. We summarize the basic results of such an analysis in 
Appendix A. Full details can be found in a companion paper [25]. 



2.5 Additional Superfields 

So far, we have introduced supersymmetry with a minimal field content of two chiral fields, X 
and <!>. Of course, there is nothing that forbids additional chiral fields, £j, from appearing in the 
action. Any sector that is decoupled from X will be protected from quadratic divergences just like 
the inflaton. It is a special feature of SUSY that we do not need to enlarge the symmetry group to 
protect additional scalar fields. We wish to treat these additional fields, = &i + V29ij)i + 9 2 Fi, 
as part of the inflationary sector. Therefore, in the rigid limit, we may write the action as 10 



Here, we assumed a decoupled action between X and Sj, just as we did for <3?. We will demand 
that (cij) = and (Fi) = during inflation. Just as in the case of <£, coupling to supergravity will 
induce masses of order H 2 for &i (unless we impose an extra global symmetry for cr,). Making 
no further assumptions, we write the action as 



U + Cx + Ux = d 4 9K(Z,^,$ + &;d^,d^,---) + /dW(£) + h.c. . (2.26) 



We have dropped the flavor index j, since there is no qualitative differences between one and 
many S's. The imaginary part of <3? is shift-symmetric, as before. 

Let us highlight some features that distinguish the action for a — eq. (2.26) — from the action 
for a— eq. (2.20): 

1) First, we note that we can have a superpotential for S, while the shift symmetry forbids 
a superpotential for <£. This implies we can write SUSY-preserving, non-derivative in- 
teractions for a. For example, the interaction A|a| 4 is obtained from the superpotential 
W = AS 3 . However, the same interaction for a is only obtained from f d 4 9 (<£ + $^) 6 
or J d 4 9 (<£ + <I>^) 4 A^A. In this case, one generates Act 4 with A oc J d 4l 9 0, where 
O = {X^X, ($ + &) 2 , • • • } is an operator whose D-term breaks SUSY. 

2) Finally, there are differences in the form of derivative couplings. For example, the coupling 
d^ipd^a can only be written by using a SUSY breaking D-term, since a and (p appear in the 
same multiplet. However, with a in a different chiral multiplet, the coupling d fl (pd^(a + a*) 
is easily achieved using /d*ef($+$t)(s_Et). This observation will be of some significance 
in our constructions of supersymmetric quasi-single-field inflation in §4.4. 

In this section, we have developed a general framework for supersymmetrizing the effective 
theory of inflation. In the next two sections, we will apply this formalism to concrete examples: in 

10 A subset of these models with £$e = was studied in [23]. For such models to leave an imprint on inflation, 
the mass of &i was required to be m|. < W~ 2 H 2 and the fluctuations in at were converted to curvature fluctuations 
£ after inflation. For a comprehensive discussion of the construction and signatures of these models we refer the 
reader to [23] . 



(2.25) 
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Section 3, we study naturalness in SUSY theories with small c s , while in Section 4, we construct 
explicit realizations of SUSY QSFI. Both sections are self-contained, so they may be read in any 
order. In particular, readers whose main interests lie in the observational signatures of SUSY are 
advised not to let themselves be distracted by the technical details of Section 3 and instead skip 
straight to Section 4. 



3 Naturalness and Small Sound Speed 

We now turn to our first application of the SUSY EFT of inflation: supersymmetric theories 
with small sound speed and large derivative interactions. As we saw in §2.1, at cubic order in 
fluctuations, there are only two single-derivative operators deforming the slow- roll action 

£ = M^H(d^f + \Ml [{d^f + l] 2 + iM| [(d^) 2 + l] 3 + • • • , (3.1) 

where 

1 2M> 4 , . 

In this section, we will realize these theories (and their generalizations to P(X) theories) in 
supersymmetry (see also [24]). In §3.1, we derive the component Lagrangians for those theories 
and show that higher-derivative couplings lead to a parametrically enhanced mass for the SUSY 
partner of the inflaton, 

H 2 

ml ~ > H 2 . (3.3) 

c s 

We then use the SUSY EFT to revisit constraints on the natural sizes of the parameters M 2 and 
M3. Specifically, in §3.2, we ask whether the hierarchy M| 3> M^l-ffj (and hence c s <C 1) is 
techincally natural. To answer the question, we construct a weakly-coupled UV-completion of 
theories with non-trivial sound speed. We then proceed to show that small c s is unnatural in 
the non-supersymmetric theory, but becomes natural when the theory is made supersymmetric. 
Finally, in §3.3, we discuss the natural size of the parameter M3 relative to the size of M2. We will 
comment on the implications of our results for the naturalness of orthogonal type non-Gaussianity, 
which relies on a cancellation between the two operators in (3.1). 

3.1 Supersymmetric P(X) Theories 

Supersymmetric sound speed. We begin with the Lagrangian for a theory with small c s , 

C Cs = -M^m 2 + H)+ M^H{d»<p) 2 + \Ml [{d^f + l] 2 . (3.4) 
In order to supersymmetrize the theory it proves useful to slightly reorganize the terms 
£c s = ( -M^H 2 + H) + \Ml ) + ( M 2 X H + \Ml ) (d^? + \ M 2 4 {d^f [{d^f + l] . (3.5) 



4 4 

2q mj 



The reason for writing the Lagrangian in this way is that (d^ip) 2 and (d^cp) 2 [(d^tp) 2 + l] each 
correspond to individual operators in SUSY, while [(d^ip) 2 + l] 2 does not. The Lagrangian (3.5) 



12 



arises from the following higher-derivative Kahler potential 
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K Cs = §(* + $t) 2 L + — d^d^] , (3.6) 

where 

c ^"l|fe = 1 + Ml-i) and C2^-^ = |(l-i). (3.7) 
The superpotential is the same as in slow-roll inflation, 

W Ca = -v{i$)X , (3.8) 

but now |^(y)| 2 = —mo- 
Curvature couplings. In our companion paper [25], we derive the complete supergravity 
action for eq. (3.6), including all couplings to the auxiliary field F, M, and b^. However, in §2.4, 
we argued that the vevs for the auxiliary fields are either suppressed or strictly zero. We therefore 
only show the result for the dominant curvature couplings [25]: 

C Cs = -\v\ 2 + Cl C l + -^-C 2 , (3.9) 

where 

A = -\d^\ 2 + (\<R-^y x \ 2 )o- 2 , (3.10) 



C 2 = - {\d^\ 2 f - 28^0-8^8^ 



+ a 2 [ ( \K - ^{Fxl 2 ) [(d^c? + (d^a) 2 ] - \U^\d^ c d^ c + d^ad v a\ 

+ \ VjNv [d^^ c d v ip c + d"ad v a] + 8^8^ + in|d M 0| 2 } . (3.11) 

Using <p c = (2M 2 i \H\) 1 / 2 t and = -2>H 2 g^ (de Sitter), we find that only four terms in (A.6) 
and (A. 7) — see underlined — contribute a mass to the field cr, 

£c s = ■ ■ ■ + [ ci (\K - H 2 ) + c 2 (±K - 2H 2 ) {d^tf + c 2 V„V v {dHd v t)\ a 2 + • • • . (3.12) 

With K = -12H 2 and V p S/ u (d ll td u t) = 9H 2 , we find 

H 2 

m 2 = QH 2 ci - \2H 2 c 2 ~ x — + 6H 2 . (3.13) 

c{ 



11 Of course, there is no unique supersymmetric generalization of these theories. There are various actions that 
differ in the couplings of fermions and additional scalars, but that reduce to the same theories when those extra 
fields have vanishing vevs. The different ways of embedding small c s -theories in supersymmetric actions were 
explored in [24]. In particular, they studied an interesting alternative operator that also generates a small sound 
speed, 

K Cs = |ci($ + $ f ) 2 + ica D^DoMa^D^' . 

As usual, the first term is only shift symmetric for the imaginary part of 4>. However, the second term in invariant 
under shifts of both the real and imaginary parts. To see this note that D Q ($ ± = D a &, using the fact that 
$ is chiral. 
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The cancellation of the leading H 2 c~ 2 terms that we find here is completely accidental and not 
protected by any symmetry. In particular, as we explain in Appendix A, in principle, there 
are a host of additional Planck-suppressed couplings between $ and X in the Kahler potential, 
cf. eq. (2.24). These terms are model-dependent, but typically they lead to large contributions 
to the mass, 

(3.14) 



ml ~ 



We therefore conclude that generically the mass of the partner of the inflaton is enhanced in the 
limit of small c s . This is important since it implies that a won't receive quantum mechanical 
fluctuations during inflation. 



SUSY P(X). The above is easily generalized to P(X) theories. Again, it is useful to first re- write the 
non-SUSY effective action 



C P{X ) = mSfoO + <{d^f + £ h<i d ^f [{d^f + 1]' 



(3.15) 



This is nothing deep, but just a simple reshuffling of the operators we are familiar with from §2.1: 



C 



P(X) 



( 4 1 4\ 



oo 



(3.16) 



-M p 2 ,(3ff 2 + ff) M^H 



Writing the effective theory in the form of eq. (3.15) has the advantage that it is easy to supersymmetrize. 
A possible form for the Kahler potential is 



K 



P(X) 



I($ + |>t)2 



l i + E h< fo*d*& + i 



n=2 



n-1 



(3.17) 



where $ = (M^\H\) 1 / 2 $. The superpotential is the same as in eq. (3.8). 



3.2 Naturalness of Small Sound Speed 

Is the hierarchy S> Mpj|ij| — and hence c s <C 1 — natural? Before we address this question, 
we digress briefly to review our prior work on the subject [40]. 

Strong coupling and new physics. It is well-known [22, 40, 41] that theories with small 
sound speed become strongly coupled at energies not too far above the Hubble scale. Of course, 
CMB observations probe modes at the Hubble scale, so, in principle, there is no problem with the 
theory becoming strongly coupled at higher energies. Nevertheless, it is interesting to take strong 
coupling as an indication for 'new physics' entering before the would-be strong coupling scale. 12 
The new physics could manifest itself as the appearance of new massive degrees of freedom or 
as a change in the physical description of the existing degrees of freedom. In [40] we applied 
this logic to inflationary models with small sound speed. We showed that the theory can stay 



12 The same philosophy can be applied to the Higgless Standard Model. A theory of massive W and Z bosons 
becomes strongly coupled around 1 TeV, unless some new physics enters before. In that case, we know that a light 
Higgs keeps the theory weakly coupled and restores perturbative unitarity. 
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weakly coupled if the dispersion relation of the inflaton field changes before the would-be strong 
coupling scale A*. We presented a specific two-field theory (see also [42-44]) that allows for large 
non-Gaussianities at Hubble without strong coupling above Hubble. However, weak coupling 
came at a price: the theory only leads to a small sound speed for unnatural mass parameters and 
couplings. 13 In this section we show that the theory, in fact, becomes natural in the presence of 
SUSY. 

A weakly- coupled UV- completion. Here, we present the basic equations of our weakly- 
coupled UV-completions of theories with small c s . For further details we refer the reader to our 
previous publication [40]. 

Consider adding a second massive field a to the slow-roll Lagrangian 



12 2 
2 m « a 



(3.18) 



where it 2 = 2M 2 l \H\ir 2 . In order for a to affect the dynamics of 7r c , we couple the two fields 
through the following interaction 



-m 3 [{d^f + l} a -> p[ir c 



1 {d^cf 



2(27^1)1/2 



(3.19) 



At high energies, u) > p, this describes two decoupled fields 7r c and a, with a small perturbative 
mixing pjr c a : 

vr a (3.20) 



At lower energies, u < p, the mixing term dominates the dynamics. The 
theory becomes a non-relativistic single-field theory (if this is counterintu- 
itive please see [40] for a detailed discussion of this feature of the theory) . 
In this regime, the disperision relation for tt is non-linear, u) = k 2 /p. For 
p 4 < M^\H\, the theory is weakly coupled at all energies up to the sym- 
metry breaking scale M 2 ^\H\. At even lower energies, uj < m^/p, the mass 
term of a becomes relevant and the theory develops a linear dispersion, 
u = c s k, with sound speed given by 



P 



(3.21) 



(K 2 il^l) 1/4 



A* 



m l/p 



If 



We see that a small sound speed requires m a <C p. Is this hierarchy 
technically natural? 

Without SUSY. The symmetries of the EFT relate the mixing term p-ka to a cubic interac- 
tion (<9 M 7r c ) 2 cr, cf. eq. (3.19). This interaction induces radiative corrections to the mass parameter 
m a and the kinetic mixing parameter p. At one loop, we find: 



13 Again, this has a parallel in the Standard Model: a light Higgs is unnatural unless something — such as SUSY- 
controls quantum corrections to the Higgs mass. 
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Sm*=* — 4 ±—<t ~ jJ^ A l ~ P 2 ' (3-22) 
dir 

5p = 7T 4 1 a ^! -p. (3.23) 

To estimate the size of the loops we have set the UV-cutoff equal to the symmetry breaking 
scale, A 4 V ~ M^\H\. We see that the weak coupling requirement, p 4 < M^ A \H\, is natural since 

2 

the loop correction in (3.23) is suppressed by the ratio 2 fgmyj < 1- However, the hierarchy 

m a < p and hence c s < 1 is not natural. Instead, the loop correction to the mass in eq. (3.22) 
implies that the natural value of the sound speed is unity, 

bc l^^^ X , (3.24) 
With SUSY. A supersymmetric version of the theory described by eqs. (3.18) and (3.19) is 



(3.25) 



The second term in (3.25) generates a quadratic kinetic mixing term with p = — ^-r . Impor- 
tantly, this breaks supersymmetry at the scale p. Notice that this SUSY breaking is associated 
with the time dependence of the background (p = t and not with the constant vacuum energy. 
Nevertheless, for p 4 < M^ A \H\, SUSY still helps to regularize the loop correction to the mass. We 
can see this in two different ways: 

1) Only loops that include the SUSY breaking will contribute to the mass of a. The loop that 
led to eq. (3.22) does not include SUSY breaking since the operator (c^7r c ) 2 cr in (3.19) is 
supersymmetric (it isn't proportional to the time-dependent background vev). SUSY will 
therefore enforce that a fermion loop cancels the boson loop leading to (3.22). Hence, we 
only get a contribution to the mass term if the SUSY breaking operator -k c a = (dot)TT c a is 
included in the loop. Such a loop then leads to 



5ml = o — \ — -a , p 2 ■ (3.26) 

dir 

We see that in the supersymmetric theory the mass of a only has a quadratic divergence, 
while the non-supersymmetric theory had a quartic divergence. This results in the maximal 

2 

radiative correction to the mass of a now being suppressed by the small ratio 2 j^. 1/2 < 1. 



16 



The hierarchy m a < p is therefore natural and protected by SUSY. 14 This implies that a 
small sound speed is technically natural 

o 2 

5c 2 ~ -A < 1 . (3.27) 

(M&\H\)V* 

We may arrive at the same result working directly in superspace by treating t f d^9 (<& + 
as a perturbation. Loops are then computed by using the superspace propagator to 
contract vertices. This leads to the following one-loop correction to eq. (3.25), 



8c loop = / d 4 e 



+ ••• , (3.28) 



where the dots indicate terms that are UV-fmite. The appearance of D 2 D 2 is the result 
of acting with the second j d 4 9 on the external lines, rather than on the propagators. 
Eq. (3.28) implies the following contribution to the mass of a, 



A 2 n 2 
A 4 p1 ' (M^Hl) 1 / 2 



5mi ~ -gM±\H\ ~ — ^— p> . (3.29) 



We again conclude that the hierarchy m a < p is natural and protected by SUSY. 

Conclusion. We have shown that the weakly-coupled UV-completion of c s <C 1 [40] is 
unnatural without SUSY, but becomes natural with SUSY. 

3.3 Naturalness of Orthogonal Non-Gaussianity 

What is the natural size of M3 ? 

We have seen that the cubic Lagrangian for P(X) theories is characterized by two indepen- 
dent operators 

L = MlMd^f + Wl [(d^) 2 + l] 2 + ^M 4 [(d^) 2 + l] 3 . (3.30) 

In the previous section, we discussed the natural value of Mi- We now turn our sights on M3. 

In [45], it was argued in the context of the strongly-coupled UV-completions, that the natural 
value of the parameter is 

, M 4 M 2 AH\ 
Ml~^~ -^J- . (3.31) 

This conclusion is of observational relevance. It implies that the interactions M|7r(5j7r) 2 and 
M^-k 3 are naturally of similar size and can therefore be cancelled against each other to produce 



14 To complete the proof, we should still show that there is no large contribution from modes with energies below 
the SUSY breaking scale p. In this regime, eq. (3.22) still applies, but the UV-cutoff is now A uv ~ p. From modes 
below the SUSY breaking scale p, we therefore get the following loop contribution to the mass of a, 

A 

r 2 P 2 

P ■ 



Ml\H\ 

This is smaller than the contribution in (3.26) and can therefore be ignored. 
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a new bispectrum shape that is orthogonal to the equilateral shape for an order one fraction of 
the natural parameter space. In the previous section, we found a natural, weakly- coupled UV- 
completion of small c s in the context of supersymmetry. We will now use this theory to determine 
the natural values of M3. In order to make contact with [45], we will determine M 4 for fixed c s 
and M^H, i.e. we will vary the parameter p, while simultaneously adjusting m a = c s p, such that 
c s remains fixed. As we discussed above, for u < c s m a , the theory describes a single degree of 
freedom, with a playing the role of the conjugate momentum of 7r, 

a^iMlmfl- 1 ^ . (3.32) 
1 m£ 

The low-energy contribution to the operator M|7r 3 therefore arises predominantly from the cou- 
pling pa 3 — i.e. 

(3-33, 

Here, M^i^l and c s are fixed, so the natural values of M3 are determined by the natural ranges 
of p and p : 

• A lower bound on p arises from the requirement that we have a linear dispersion, uj = c s k, 
at horizon crossing. This means that c s m a = c 2 s p > H. An upper bound on p arises from 
the lower bound on m a , associated with the minimal loop contribution in eq. (3.26): namely 
ml > 5ml ~ p 2 (M^\H\)~ 1 / 2 p 2 leads to p < c(Mjlfll) 1 / 4 . Therefore, we find 



-2 



rrH < p < ciM^lHl) 1 /* . (3.34) 



The natural range of p can be determined by simple loop estimates, as before. The 
lower bound is determined from the loop contribution to p arising from the interaction 
p(M 2 l \H\)~ 1 / 2 (d /1 7T c ) 2 a. As before, the coupling only gets renormalized if we insert the 
SUSY breaking operator pn c a inside the loop 



** — — ^ffi ■ < 3 - 35 » 

Naturalness of the coupling then requires p > 5p ~ p 5 (Mp 1 |ij|) _1 . An upper bound on p 
follows from the renormalization of the mass of a, i.e. 5ml ~ p 2 . Because we have fixed 
m a in terms of c s and p, we must have p < m a . Therefore, we find, 



P 5 < < 



^p 2 il# 

Combining the ranges of p and p with (3.33), we find that the natural range of M| is 



(3.36) 



—Rll — !. < Ml < Hi 



< Ml < ±- , (3.37) 
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where A c = J fiT 2 /(4M pl |i/|) 1 / 2 ~ 10~ 5 . We see that we can naturally make the value of M3 larger 
than the natural value in the strongly-coupled theory [45], but not smaller. 

Conclusion. We have shown that the natural range of M3 is larger in a weakly-coupled 
SUSY completion than in strongly-coupled models. Interestingly, the lower bound on M3 is 
identical in both cases, Mg ~ M^HlcJ 4 , but the upper bound is now higher. The orthogonal 
shape was shown in [45] to arise for an order one range of parameters in the strongly-coupled 
theory. The fraction of the natural parameter space that gives orthogonal shape is likely to 
decrease in the weakly-coupled SUSY model. 

4 Supersymmetric Quasi-Single-Field Inflation 

One of the most exciting applications of the supersymmetric effective theory of inflation is to the 
quasi-single-field inflation (QSFI) models of Chen and Wang [11]. These models involve extra 
fields besides the inflaton field, with masses close to the Hubble scale. This feature makes su- 
persymmetry a natural arena for QSFI. In this section, we present a supersymmetric completion 
of the model of [11]. Moreover, we use the effective theory approach to propose some obvious 
generalizations of the original model. We then determine which of these new possibilities can nat- 
urally arise in SUSY. As advertised in the Introduction, QSFI allows for large non-Gaussianities 
with a unique signature in the squeezed limit. We will show that this signature is robust and not 
sensitive to details. 

4.1 Quasi-Single-Field Inflation 

The basic idea behind QSFI is very intuitive. In this section, we will explain the mechanism and 
its key predictions in simple physical terms. We then provide a systematic classification of all 
possible variations of QSFI and show which of them are naturally realized in SUSY. 

Basic Mechanism and Perturbative Predictions 

Setup. QSFI couples a massless Goldstone mode, £ 7r = — ^(d^itc) 2 , to a 
massive isocurvaton mode, C a = — ^(d^a) 2 — ^m 2 a 2 , through a quadratic 
mixing term £ m ix, which [11] chose as £ m ix = p^cC. Incidentally, this is the 
same mixing term that featured prominently in our weakly-coupled UV- 
completion of theories with small sound speed, cf. §3.2. This coupling arises 
if the background fields follow a curved trajectory with constant radius of 
curvature and constant angular velocity. The mixing converts fluctuations in a into fluctuations 
in 7r c and hence (. If the mass of the second field is at or below the Hubble scale, m a < 
then quantum fluctuations in a can contribute significantly to the final curvature perturbation. 
Moreover, interactions in the cr-sector are much less constrained than interactions in the inflaton 
direction. A large interaction term £; nt in the cr-Lagrangian can then be an important source of 
non-Gaussianity. The specific example £j n t = — /xcr 3 — Act 4 was explored in [11]. 

Power spectrum. For p < H, we can treat the mixing term as a perturbative correction, 
while the leading order dynamics is determined by 

Co = -\{d^ c ) 2 ~ |(<V) 2 - \mla 2 . (4.1) 
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The quadratic mixing is then captured by the following transfer vertex 



C 



mix = P^c<y 



7T 



-0- a . 



(4.2) 



The leading perturbative correction to the power spectrum of £ is 

1 H 4 



H0 0- 



O A? = - — — 
4 4Mj|#| 



l + c W (| 



(4.3) 



where the precise form of the function c{v) can be found in [11]. Since the correction is scale 
invariant, it only leads to an unobservable shift in the overall normalization of the power spectrum. 



Scale invariance. It is worth digressing briefly to explain why QSFI does not introduce large 
violations of scale invariance. Since the mixing term couples the Goldstone boson 7r to a massive field a, 
one might suspect that the Goldstone boson effectively becomes massive, leading to a violation of scale 
invariance in the power spectrum of order p 2 /H 2 . Moreover, one may worry that the violation of scale 
invariance in the c-sector, of order rri^/H 2 , gets communicated to the 7r-sector. However, the result of the 
explicit calculation in (4.3) shows no violation of scale invariance. What happened? 

Exact scale invariance of the power spectrum is the result of an exact global symmetry under which 
t — > t + c and k — > ke Hc (in de Sitter space). The transformation on t can be undone by a time 
diffeomorphism, after which one is left with an equivalent global symmetry tt — > tt + c. If this symmetry is 
unbroken, then every mode experiences the same history and the power spectrum of tt is scale invariant. 15 
In QSFI, the action for tt and a posses such a symmetry and hence there is no reason to expect any 
violation of scale invariance. 

This explanation, while true, does not fully address the concern regarding massive fields. If the 
isocurvaton fluctuations, a, were converted into curvature perturbations, £ = — Htt, after inflation ends, 
then we would indeed find violations of scale invariance. This is because the symmetry t — > t + c is 
also broken in order to end inflation. The evolution of a mode that is frozen outside the horizon is 
insensitive to the end of inflation and no violation of scale invariance appears. However, for modes that 
do not freeze out (i.e. massive modes), the conversion to curvature perturbations at the end of inflation 
picks out a scale, since the amplitudes of the different modes are measured at a specific time. For QSFI, 
the end of inflation does not induce a significant violation of scale invariance, because the isocurvature 
perturbations are converted to curvature perturbations before inflation ends. The curvature perturbations 
£ are constant outside the horizon. By converting a into £ during inflation, the observable £ modes are 
therefore insensitive to the end of inflation. 

Bispectrum. The most intriguing aspect of the phenomenology of QSFI is the fact that 
it leads to large non-Gaussianities with a unique scaling behavior in the squeezed limit. As we 
alluded to in the Introduction, this provides the opportunities to use CMB and LSS measurements 
to probe Hubble-mass degrees of freedom during inflation. 

Chen and Wang [11] explicitly computed the three-point correlation function induced by 
the interaction £ mt = — ^c 3 . In Appendix B, we show how to reproduce the most important 
features of their answer from simple physical reasonings and back-of-the-envelope estimates. For 

15 In most modeis this globai time transiation symmetry is oniy approximate, being weakiy broken by the time 
cvoiution of the Hubble scale, H ^ 0. 
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the amplitude of the bispectrum we find 



Ant = o & / NL ~ A" 1 (|) 3 1 . (4.4) 



Notice the enhancement of /nl by the inverse of the amplitude of the primordial fluctuations, 
A^ 1 ~ 10 5 . This allows for a large non-Gaussianity even in the perturbative regime with p < H 
and p < H. In Appendix B, we also state simple "Feynman rules" that allow us to estimate the 
size of more general diagrams. Essentially, each mixing insertion contributes a factor of p/H and 
each interaction vertex gives a factor of p/H. The powers of A^ 1 are determined by the number 
of a's in the interaction vertex, with interactions containing fewer <r's being suppressed. 

The a field is massive and decays outside of the horizon as ( — t)*^ u , where v = \J | — -jp 
and r is conformal time. This leads to a non-trivial squeezed limit for the bispectrum 

hm (C kl Ck 2 Ck 3 ) oc • (4.5) 

This scaling behavior is intermediate between that of the local shape (k^ 3 ) and that of the 
equilateral shape (k^ 1 ). By measuring the squeezed limit we determine the index v and hence 
the mass of the isocurvaton a. 

The result (4.5) has a simple physical interpretation (see Appendix B): First, we recall that 
the squeezed limit corresponds to the correlation between a long-wavelength mode k\ and two 
short- wavelength modes &2 ~ The long mode crosses the horizon at \k\T\\ ~ 1, some time 
before the horizon crossing of the short modes at \k2Ti\ ~ 1. The superhorizon evolution of the 
long mode between t\ and T2 leads to a suppression of its amplitude 

0*1 fa) ~ cr^fa) ( ^rj ~ o-feifa) ( y 2 ) ■ (4.6) 



We can write eq. (4.5) in the following suggestive way, 



3/2-1/ 



1 { k \' 

^ k ^ )K ||U) • (4 - 7) 

local 

We recognize this as the product of the local shape and a modulation that precisely matches 
the suppression of the long mode in eq. (4.6). This makes a lot of sense: if the isocurvaton is 
massless it freezes after horizon crossing and we expect local type non-Gaussianity from the non- 
derivative interaction a 3 . The squeezed limit of QSFI can therefore be interpreted as a modulated 
local shape with the modulation determined by the superhorizon evolution and hence the mass 
of the isocurvaton mode. 

Future data [15, 16] will have much to say about the primordial bispectrum in the squeezed 
limit. Here, we want to highlight the role of LSS observations. In recent years, the scale- dependent 
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bias has emerged as a sensitive probe of primordial non-Gaussianity [17]. Receiving most of its 
signal from the squeezed limit of the bispectrum, the scale-dependent bias is an ideal probe of 
QSFI. Eq. (4.5) implies the following scaling for the non-Gaussian bias 

This prediction motivates generalizing the LSS data analysis to include v as a free parameter; 
whereas, to date, the analysis has mostly been restricted to the case Ab oc k~ 2 [19] (but see [13, 
20]). 

Trispectrum. QSFI makes further interesting predictions for the four-point function: First, 
we see that, correlated with the three-point function from the a 3 interaction there is a four-point 
function from a scalar exchange diagram 



Ant = -^ 3 > y — A & r NL ~A- 2 (-|) (it) »/ 2 L . (4.9) 

/ \ 

We note that QSFI gives a natural mechanism to boost the amplitude tnl relative to (|/nl) 2 - 
(Recall that if a single source is responsible both for the power spectrum and the non-Gaussianity 
of ( then tnl = (|/nl) 2 -) We expect this feature to lead to a stochastic halo bias (i.e. bias 
inferred from (6^5^) ^ bias inferred from (<5h<W) [46, 47]. These signatures of QSFI deserve 
further investigation [48]. 

For models with additional quartic couplings such as cr 4 , we get a four-point function asso- 
ciated with a contact interaction 



C U!l --Xa i X <?NL~ A" 2 (|:) 4 A . (4.10) 



The shapes of the four-point functions in (4.9) and (4.10) are distinct. The scalar exchange non- 
Gaussianity peaks in the collapsed configuration where all momenta have equal amplitude, while 
the contact interaction peaks in the squeezed configuration where one of the momenta vanishes. 

We will return to the phenomenological predictions of QSFI in future work [48]. In this 
paper, we will instead explore the theoretical foundations of QSFI and its relation to supersym- 
metry. 

Variants of Quasi-Single-Field Inflation 

The most general form of quasi-single-field inflation can be written as a deformation of slow-roll 
inflation: £qsfi = £ s . r . + A£, where 

AC = -^(<V) 2 - \mla 2 + £ mix + £ int . (4.11) 
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It is straightforward to extend the ideas of [11] to different mixing terms £ m i x and/or different 
interactions £; n t- In this section, we will be systematic about this and discuss the range of 
possibilities in the effective theory approach. We begin with a non-SUSY treatment and then 
discuss its SUSY implementation. 



Mixing terms. 
mixing terms 



Instead of C $ = m^Kd^ip) 2 + l]a — > ptt c cj we could consider the following 



,(2) 
'mix 

(3) 



m 



a,/3 



/{M'^HD 1 / 2 . However, it is easy to see that £ mix does not lead to observable 



air c a , 

(3) 



(4.12) 
(4.13) 



where {a, (3} 

signatures. Integrating by parts, we find £^ lx — > (V 2 7r)<7 + V^(V At 7rcr) — 3Ha]. In the 

decoupling limit, H — > 0, we have V 2 7r = and therefore, the mixing term will not contribute. 
Equivalently, we may remove the mixing term by a field redefinition tt c — > ir c — Because 
a — > as t — > 0, such a field redefinition does not change the late-time correlation functions. For 
these reasons, we may set j3 = without loss of generality. 16 

(2) 

Introducing £ mix will modify the dispersion relations of physical modes. We will require 
that no group velocity exceeds the speed of light. Consider the Lagrangian 



£kin = \ (tt c cr) 
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(4.14) 



where we have allowed for sound speeds for both 7r and a: i.e. c s 2 = 1 + and c s 2 = 1 + e CT . 
Solving the equations of motion, we find two positive frequency solutions with velocities given by 



1 



v± 



1 + Ue a + e7r ) ± 



) 2 + a 2 



(4.15) 



The constraint v + < 1 requires that e a e n > a 2 . If we take e CT = e n = e, the velocities simplify to 
v± = (1 + e ± a) -1 , which makes e > a transparent. In the perturbative regime, a < 1, we only 
require a small deviation from the speed of light, c = 1. The fact that we have to introduce non- 
trivial sound speeds for both n and a to avoid superluminal propagation is completely consistent 
with expectations from quantum field theory in flat space. Note that sound speeds for ir and a 
arise from operators which include terms of the form (d^d^ir) 2 , (d^ad^n) 2 , etc. It was shown 
in [49], that these types of four-derivative operators have coefficients that are constrained by 
unitarity to have non-zero values. Although in the effective theory, it appears that we have to 
adjust these coefficients by hand, any manifestly Lorentz-invariant UV-completion of the effective 
theory will necessarily yield the coefficients consistent with the constraints. 



In summary, we will consider two mixing terms £ mix = p~K c <y and £ mix = a-h c a. For the 
mixing to be a perturbative effect, we require {jj, a} < 1. Note that this is not a model-building 
requirement, and, in principle, we can allow the mixing parameters to be order one or larger [48]. 



We thank Kendrick Smith for helpful discussions on these issues. 
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Interactions. These two quadratic mixing terms may be combined with any of the following 
cubic interactions 





= m 3 a [(d^) 2 + l]a 
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(4.16) 
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(4.27) 


£-8b 
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£-9a 


= A 3 2 a(c» 2 


(4.29) 


£-9b 


= Af(d^a) 3 
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£-9b 


= A 4 " V 


(4.30) 



This list captures all possible combinations of the fields ip and a. We have restricted to operators 
with at most one derivative acting on each field, but one may easily include additional derivatives. 
In Appendix B, we show how to estimate the size of non-Gaussianities for these interactions. The 
results are summarized in Table 1. In the table we also indicate which interactions, in principle, 
allow for large non-Gaussianities (large NG), which peak in the squeezed limit (S.L.), which 
preserve supersymmetry (SUSY), and which ultimately lead to natural models (Natural). 

In the remainder of this section we will explore which of these new models of QSFI have 
a natural microphysical implementation. We will distinguish cases that are strictly unnatural 
because they involve fine-tuning at tree level (§4.2); cases whose naturalness remains to be estab- 
lished because they are likely to involve fine-tuning at loop level (§4.3); and finally cases that are 
completely natural because supersymmetry regulates the radiative corrections (§4.4). 

4.2 Problems at Tree Level 

When organizing the effective theory of inflation in §2.1, we found it convenient to expand the 
action in terms of operators like [(d^cp) 2 + 1]. However, when discussing the natural values of 
parameters, one should think of (d^(p) 2 and 1 as independent operators. In general, writing 
terms like [(d^ip) 2 + \\0((p), without independently adding the operator O(p) should be viewed 
as a tree-level fine-tuning. 17 This is most apparent in the context of supersymmetry, where the 

17 The fine-tuning of terms necessary to cancel tadpoles is an exception to this rule about naturalness. In 
writing the EFT of inflation, we demanded that our background solution takes a specific form and then fix certain 
coefficients to cancel tadpoles. Tadpole cancelation is not a tuning of the action, but simply represents the 



24 



Table 1: Non-Gaussianity in Quasi-Single Field Inflation. 
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operators proportional to 1 and (d^ip) 2 typically arise from fundamentally different terms in the 
Kahler or superpotential, cf. §3.1. Cancelling the constant contribution from (d^<p) 2 only occurs 
when the coefficients of a priori independent operators are carefully chosen to cancel. With this 
in mind we wrote the operators in eqs. (4.16)-(4.30) without the combination [(d^ip) 2 + 1], unless 
it was related to tadpole cancellation in either ir or a. 

It then follows immediately that the interactions L\ - £5 can't lead to large non-Gaussianities, 
unless the tree-level action is fine-tuned. Let us demonstrate this case-by-case: 

• The interactions £i a ,& come from the same operator £1 that leads to the mixing term 
pTT c °~-, where p = w? a b / {M 2 ^\H\) 1 / 2 . We therefore can't make the interaction large without 
inducing a large mixing term. This constrains the non-Gaussianity to be small, /nl < 1- 
For £ift, one may decouple the interaction and mixing terms, but this requires fine tuning. 

freedom to choose initial conditions. However, beyond the cancelation of tadpoles, the cancelation of the constant 
in (d^ip) 2 = — 1 + ■ ■ • will typically require tree-level fine-tuning. In some exceptional cases, there may be a 
dynamical explanation for the cancelations. However, in general, we should be suspicious of exact cancelations 
between operators. 
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• Similarly, the interaction £2 comes from the same operator that leads to the mixing term 
ocrr c &. Again, this constrains the non-Gaussianity to be small, /nl < 1- 

• The interaction £3 could give large non-Gaussianity for rh 3> H. However, the operator 
£3 = rh 2 {d ^ip) 2 a 2 C m 2 a 2 also contains a tree-level mass term for a. Without fine-tuning, 
we require rh 2 < H 2 and hence /nl < 1- 

• Similarly, the interactions £4 a ,& are necessarily related to the operator fh a> i,aa — > Hfha^a 2 . 
Keeping the mass naturally small requires m a ^ < H and forces the non-Gaussianity to be 
small, /nl < 1. 

• The interactions C^ a ^ )C only lead to large non-Gaussianities if \ a ,b,c S> 1. Again, it is 
straightforward to see that this requires fine-tuning. The interactions C^ a ^^ c come from 
operators that also contain a contribution to the kinetic term of the form A a b,c& 2 - Without 
fine tuning, if \ a ,b,c 3> 1 the true canonical field is a c ~ c/A 1 / 2 and /nl "C 1. 

On these grounds we reject the interactions £1 - £5 as candidates for natural models of QSFI. 
For the remaining interactions naturalness is a bit less trivial to check. 

4.3 Challenges at Loop Level 

Naturalness criteria. Given an EFT with a set of couplings, one should ask if these parameters 
are stable under radiative corrections. Loop contributions to various parameters can be estimated 
as a function of the UV-cutoff of the theory. If a parameter receives large corrections, we call the 
parameter fine tuned. In many scenarios, the UV-cutoff of the EFT is unknown and, in principle, 
may be taken to much higher energies than are probed directly in experiments. Parameters may 
become unnatural if the UV-cutoff (and hence the range of validity of the EFT) is taken to be too 
large. One may also reverse the logic to predict the breakdown of the EFT based on naturalness. 

Before discussing supersymmetry, it will be useful to first understand the naturalness of 
QSFI with minimal field content and interactions. Models based on the interactions C% and 
£7 are effectively described by a slow-roll background with weakly-coupled interactions, so, in 
principle, no new physics is required up to the symmetry breaking scale Aj, ~ (M 2 [|£T|) 1//4 [40]. 
In that case, the UV-cutoff of the effective theory can be as large as A uv ~ Aj>. On the other 
hand, models based on £$ and £9 may become strongly coupled at energies below the symmetry 
breaking scale. As we discussed in §3.2, this can be viewed as an indication that new physics 
should appear at or below that strong coupling scale, A* = Aj. Hence, for models based on £g 
and £9 we impose A uv = min{Ab, Aj}. 

The absence of any additional new physics up to Aj, or A* leads to a relatively strong 
condition for naturalness: namely, the theory ought to be natural for A uv ~ min{A{,, A*} 3> H. 
However, of course, we can't exclude the possibility that new physics may appear below the scale 
min{Ah, A*} and that this would improve the radiative stability of the models. Ultimately, SUSY 
will play this role for some of the interactions. 
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Loop corrections. A principle worry is that the large interactions required for observable 
non-Gaussianities induce large loop corrections to the mass of the isocurvaton: 

5ml = \ \ . (4.31) 

\ / 

•v. S 

For the QSFI mechanism to be operative, we require m\ < \H 2 . A necessary condition for 
naturalness is therefore 5ml ~ \H 2 . Let us study the interactions C§ - Cg one- by-one: 

• The interaction is somewhat special, in the sense that loop contributions to the mass 
are finite 

Sm 2 a{6} ~ /x 2 . (4.32) 

This is smaller than the supergravity-induced mass ml ~ H 2 , if \i < H. Hence, Cq is a 
promising candidate for supersymmetric QSFI. Indeed, in §4.4, we will present a natural 
model based on Cq. 

• Without SUSY, the interactions £7 - Cg all give UV-divergent loop corrections: 

- The interaction £7 induces the following one-loop mass renormalization 

5ml {7} ~ A 2 A 2 V . (4.33) 

In the absence of additional new physics below the symmetry breaking scale, we cut off 
the loop at A uv ~ A^ ~ ( Af 2 [ I H\ ) Requiring the loop contribution to be subdominant 
relative to the tree-level contribution, m a ~ H, puts a bound on the size of the non- 
Gaussianity 

/ NL{7} < a 3 A- 1/2 . (4.34) 

This allows marginally observable non-Gaussianity, but only under the optimistic assump- 
tion a ~ 1. 

- The interactions Cg leads to the following radiative contribution to the mass 

1Y 1,2 

where Ax,2 > H in order for the theory to be weakly coupled at horizon exit. The naturalness 
condition 5m a ^ 8 y < H, therefore requires A uv < Ai,2. This cutoff is only consistent with 
A uv = min{A(,, Ai^} for Ai : 2 > A& and A uv ~ A&. Naturalness then strongly constrains the 
non-Gaussianity associated with Cg, 

/nl{8} < « 3 . (4.36) 

Hence, without a way to regularize the loops below the symmetry breaking scale, the 
interactions C% a ^ don't naturally lead to observable non-Gaussianities. Alternatively, we 
can keep a more open mind about the possibility of new physics appearing at relatively low 
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energies. To achieve a given /nl with natural parameters, one then finds a bound on the 
UV-cutoff 



a 



,3/2 



• {m 2 Jh\) 1/a . 




1/2 



(4.37) 



fi 



NL 



- The interactions Cg leads to wavefunction renormalization 



iv 3,4 



(4.38) 



For A uv ~ minjA^, A3.4}, this is always smaller than the tree-level value Z a = 1. The 
interactions £9 can therefore lead to natural models with large /nl- However, being pure 
derivative interactions, the non-Gaussianity is suppressed in the squeezed limit, making 
these models less interesting for our present purposes. 

4.4 Supersymmetry and Naturalness 

Supersymmetry plays two important roles in the microphysical implementations of quasi-single- 
field inflation: 

1. SUSY naturally and inevitably introduces at least one additional scalar field with mass of 
order the Hubble scale H. 

2. SUSY can help to regulate dangerous loop corrections that would otherwise become impor- 
tant in the limit of large interactions and threaten the naturalness of the model. 

Can SUSY help to regulate the UV- divergences that we found in the previous section? Only 
for £§ a . It is easy to understand why: as we discussed in §2.2, spontaneous breaking of Lorentz 
invariance necessarily breaks supersymmetry. Therefore, any interaction that is not manifestly 
Lorentz invariant must be proportional to SUSY breaking. This is the case for the cubic interac- 
tions £7, £$b and Cg a ,b- Therefore, when we embed these interactions in a supersymmetric theory, 
there will be no extra cancelations below the SUSY breaking scale, namely u> < = (M^ l \H\) 1 ^' 1 . 
What distinguishes C% a is that this interaction does not know about SUSY breaking directly. The 
same interaction can occur in theories with unbroken SUSY. Because the mass for a arises from 
SUSY breaking, it will be further suppressed relative to the estimate in (4.35). Although the 
interaction £g a preserves SUSY, the mixing terms pir c a and cefr c & don't. A contribution to the 
mass of a is therefore generated by inserting the mixing terms inside the loop. 18 The associated 
loops are still divergent, but suppressed relative to (4.35). For the mixing a7r c cr, the suppression 
is rather mild, since the order of the divergence stays the same, 




5ml 



a{8a,a} 




f a ~ 



uv 



(4.39) 



18 



This is the same effect that we discussed in §3.2; see eq. (3.26). 
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In the regime of interest, Ai < A& and hence A uv ~ A%. Naturalness, Sm^^^y ~ aAi < H, then 
puts a bound on the size of non-Gaussianity 

/NL{8a,a} < A^ 1 . (4.40) 

This could marginally be observable, but only in the questionable limit that Ai get dangerously 
close to H. 

It is more promising to work with the mixing term p7r c cx. In this case, the order of the 
divergence is reduced, 

Sm l{8a,p} = V---4 9 ^~ T2 A uv ( 4 - 41 ) 

^ 

This suppression is sufficient to make the model natural. In particular, for A uv ~ min{A{,,Ai}, 
we find 5rn a {$ a } < p < H. Naturalness therefore doesn't put an additional constraint on the size 
of non-Gaussianity. In particular we can get large non-Gaussianity even for Ai safely above H, 

/NL{8a,p} < ^A- 1 . (4.42) 

— 1/2 

Even for Ai — > A&, the non-Gaussianity can be observably large, /NL{8a,p} < A^ . The interac- 
tion £$ a is therefore a promising candidate for a natural SUSY model of QSFI. 

Although SUSY does not regulate the divergences associated with £7 and £gfc; there is 
no direct constraint on /nl if the cutoff is lowered. We may have still have natural models 
for these interactions, but it would require new physics to enter at low energies. We distinguish 
these situations from 'natural SUSY implementations of QSFI', as they require something beyond 
SUSY (or in place of SUSY) to be natural. 

We have identified two interactions £i n t = {Cs^C&a} with the potential to lead to natural 
models under minimal assumptions. Next, we will discuss whether these models can indeed be 
consistently embedded in a supersymmetric framework. 

Model I 

We first present a SUSY implementation of QSFI using the interaction Cq. In the interest of 
economy, we would like to find a model where the massive field, a, is part of the same chiral 
multiplet as (p. This choice restricts the form that the model can take. We can't obtain Cq from 
a superpotential, since that would break the shift symmetry ip — > tp 4- const. Derivative mixings, 
such as d^TTcd^a, are also difficult to generate, making it most natural to consider the mixing 
term p7r c a. 

Making these choices, the SUSY implementation of the model with £ m i x = p7r c °~ and £j nt = 
Cq = —ua 3 is 

4!sfi = Jd i e[l($ + &) 2 + ^ + &f + ^(<S> + &) 5 ] +C X , (4.43) 

' V ' > V ' 
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with 



P = \ — and p = f ^\ 3 ' . (4.44) 



Ax -~ (A 



The constraints p < H and p < H, require Ai > ^M p \ and A 2 > e 1//3 (^) 1/3 M p i ~ 
(A^) 1 / 3 y / eM p i. We see that Planck-suppressed corrections are sufficient to generate the re- 
quired couplings. An alternative way to generate the cubic coupling is the operator K in t = 
tAj($ + &fX^X. In this case p ~ M 2 l H 2 /(A 2 ) 3 , which implies /i < # if A 2 > {eA^^M^. 
We have not written an explicit mass term for a since it will be generated by supergravity cor- 
rections: m a ~ H. Supergravity will also generate a cubic term, but with a coefficient that is 
suppressed by Ai, namely p ~ 2>H 2 / k\ <C H. 

As we noted above, even in the absence of SUSY the mass of a is does not receive dangerously 
large corrections, if p < H and p < H. Embedding the model in supersymmetry does not change 
this situation. Explicit loop calculations in the supersymmetric model simply reproduce the 
results of the previous section. The role of SUSY in this model is to provide an explanation for 
the extra Hubble-mass scalar a. 

Model II 

Next, we discuss a natural SUSY implementation of our second candidate model, Cg a . As in the 
previous model, it would be ideal to embed a in a single chiral multiplet with tp. We will start 
with this assumption, only to find it cannot lead to large /nl- However, we will show that adding 
a second chiral field, X, leads to natural models with large non-Gaussianity. 

In order to generate /nl > 1 from the interaction Cs a , we require 

Ai < HA^ 1 ~ A-^iM^Hl) 1 / 4 . (4.45) 

In this case, SUSY is required to regulate the loop correction to the mass of a. Indeed, we can 
write the operator £$ a supersymmetrically 

j d A 6^-(§ + $t)3 c —aid^a) 2 . (4.46) 

Since SUSY is unbroken by the interaction, the previously dangerous loop can now be controlled. 
However, eq. (4.46) necessarily includes the mixing term 

M 2 H 

£mix = Hr— 7T<r • (4.47) 

With the constraint (4.45), this leads to p^cC, with p S> H. This spoils the QSFI mechanism, 
since large p implies a large effective mass for a. We therefore find that QSFI with interaction 
£g a can't be realized in a minimal setup with only a single chiral superfield 

It is clear that this obstruction to naturalness is specific to the case of a single chiral field. 
Next, we explore the possibility of realizing the model with additional chiral fields. It is now 
straightforward to separate the mixing term pir c a from the interaction C% a . It is easy to see that 
the following Lagrangian leads to the desired couplings, 

£g } FI = j d 4 0[|($ + $t)2 + EE t + l.($ + $t)2 (E + s t) + _L( E + S t) ES tj (4.4 8 ) 

S v ' ^ v ' 

^-mix -^int 
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This action contains a tadpole for a that we are implicitly canceling through an additional small 
coupling to X. The spectator field S does not break SUSY by itself, but only indirectly through 
its coupling to the inflaton <I>. In the limit A — > oo, the action for X decouples and describes 
a supersymmetric subsector (prior to coupling to supergravity, as usual). By the usual SUSY 
non-renormalization theorems, we will not generate a mass for a through any loop diagram. For 
example, if we were to insert (K- mt ) 2 and compute the loop, we would indeed find a divergence. 
However, the divergence is only logarithmic and contributes to the wavefunction renormalization 
Z f d^OTiY^ rather than to the mass term. By inserting (-ftT m i x ) 2 , we can introduce a mass for a, 
but this mass is smaller than Hubble when we cut off the loop at A uv ~ min{Ab, Ai}. The theory 
described by eq. (4.48) is therefore technically natural. 

4.5 Summary 

We summarize our conclusions regarding the naturalness of quasi-single-field inflation in the 
presence of supersymmetry: 

• C\ - £5 are unnatural and require fine-tuning at tree-level. 

• £7 and Csb require additional physics beyond SUSY for naturalness. 

• £g a ,6 are natural, but the signal is suppressed in the squeezed limit. 

• Cq is natural. Loops are finite even in the absence of SUSY, but SUSY explains the origin 
of the massive isocurvaton. The model requires only a single chiral multiplet that contains 
both the inflaton and the isocurvaton. 

• £g a is natural, but only if the inflaton and the isocurvaton live in separate multiplets. SUSY 
is required to explain radiative stability. 

5 Conclusions 

All of the current CMB and LSS data is completely consistent with the simplest effective field 
theory of inflation — a theory of Goldstone bosons of spontaneously broken time translations [22]. 
However, despite this phenomenological success, the fundamental origin of the inflationary ex- 
pansion remains a mystery. The basic challenge is easy to understand: it is hard to protect 
fundamental scalar fields from radiative corrections. The simplest field theory models of inflation 
are therefore unnatural in the technical sense [3]. Taking this naturalness problem seriously, one 
is led to consider symmetries as a way of protecting the inflationary background from destructive 
quantum corrections. This means either internal global symmetries or supersymmetry. The role 
of global symmetries in inflationary models is well-understood, although generic Planck-scale 
breaking effects are rarely included [36, 37]. Unfortunately, it seems that models with global 
symmetries are observationally indistinguishable from models that are simply fine-tuned (but 
see [23, 50]). In the case of supersymmetry, on the other hand, the observational prospects are 
considerably more optimistic. Supersymmetry requires that the (real) inflaton field has a scalar 
partner in order to match the two degrees of freedom of the fermionic superpartner. When SUSY 
is spontaneously broken by the spacetime curvature of the de Sitter background, this second 
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scalar receives a mass of order H (while the inflaton stays light either due to an additional global 
symmetry or due to some moderate fine-tuning). This feature is special to SUSY, as internal sym- 
metries cannot be spontaneously broken by curvature without violating the Coleman-Mandula 
theorem [51]. We have proposed to use this massive scalar as a window into the supersymmetric 
origin of inflation. Fortunately, the presence of a Hubble-mass scalar during inflation has a clean 
observational signature in the squeezed limit of the primordial bispectrum [11]. Notice that the 
path from naturalness via SUSY to Hubble-mass fields and the squeezed bispectrum requires 
little more than simply following your nose. 

An important by-product of this paper was the development of a general framework for 
systematically incorporating supersymmetry into the effective theory of inflationary fluctuations. 
This supersymmetric effective theory of inflation has a wide range of applicability. We presented 
two examples in detail: models with small sound speed and slow-roll quasi-single-field inflation. 
We showed how supersymmetry helps to address naturalness issues in a faithful way. It would be 
interesting to study the broader phenomenology of the SUSY EFT of inflation: models with higher 
derivatives, quasi-single-field inflation with small sound speed, etc. Moreover, our treatment 
made some simplifying assumptions concerning the effects of SUSY breaking on the inflationary 
dynamics (see §2.2; in particular, we chose separate multiplets for the SUSY breaking spurion X 
and the inflaton <£.). We also restricted the size of superpotential terms and hence the couplings 
to the auxiliary supergravity fields. Finally, we didn't study scenarios in which the fluctuations 
in additional isocurvature fields are converted to curvature perturbations after inflation (see 
e.g. [23]). Relaxing those assumptions will allow a broader class of models. This may lead to 
additional observational signatures that haven't been anticipated yet. In the meantime, we should 
also determine how present data constrains the models presented in this paper. Efforts are under 
way to see what CMB and LSS data has to say about SUSY in the sky [48]. 
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A Supergravity for Effective Theories 



In a companion paper [25], we study how general effective theories with global supersymmetry 
are coupled to minimal supergravity. In this appendix, we give a brief, and mostly qualitative, 
overview of that work. 

The basic challenge in constructing a theory of supergravity arises from the fact that rigid 
SUSY is a spacetime symmetry. Specifically, the SUSY algebra {Q a ,Qa} = 2(7^ P M relates 
SUSY transformations to spacetime translations. In order to describe a supersymmetric theory 
on a curved background, one needs to promote both vector and spinor derivatives (c^, D a , and 
D a ,) to covariant derivatives (V^, V a and V a ). To do so, one promotes the metric and the spin 
connection to superfields, and imposes constraints and gauge-fixing conditions that reduce the 
number of new, off-shell, degrees of freedom to a manageable number. Minimal supergravity 
represents one set of self-consistent constraints that allow us to define a covariant theory [39]. 
The off-shell description contains two dynamical fields, the vierbein and the gravitino ip", and 
two auxiliary fields, a complex scalar M and a real vector b^. All geometric objects may be 
defined in terms of these fields. Superspace is then defined with new fermionic coordinates O. 
With some work, the action in superspace may be written as 



- 2 J d 4 x| d 2 eiJ§(P 2 -8i?)e-T^ X ' xW ^^ 



+ h.c. , (A.l) 



where k = 1/M P \. There are two crucial elements that make this action manifestly supersymmet- 
ric: First, the operator (T> 2 — 8R), with curvature superfield R D —\M + j^IZQ 2 , projects any 
scalar operator O onto a chiral operator C = (T> 2 — 8R)0 (i.e an operator satisfying T> a C = 0). 
Second, the chiral density, E D det e (1 — 2 M) , is defined such that J d 2 @ E C is supersymmetric 
for any chiral operator C. 

In inflation, the most important supergravity contributions are those that arise from cur- 
vature couplings. In minimal supergravity, the coefficients of the curvature couplings cannot be 
changed arbitrarily and they will induce masses and couplings that are controlled by the Hubble 
scale H. For example, for slow-roll inflation, with A" s .r. = + + XX* , eq. (A.l) leads to 
curvature couplings of the form 

£su g ra3-^exp[-^(a 2 + |x| 2 )] , (A.2) 

where 1Z = —12H 2 is the curvature of the de Sitter background. When K(<&, X, X\ . . .) = 
K($, $t, X, At) = K($i, $J) (i.e. when there are no higher-derivative terms or additional curva- 
ture couplings), the curvature coupling arises from a universal coupling to the chiral curvature 
superfield R. In this case, one normally Weyl rescales the metric to go to Einstein frame and 
these curvature terms then appear as part of the supergravity scalar potential 

mWra = e* K fo>to (g'WiWDjW - 3k 2 \W\ 2 ) , (A.3) 

where DiW = d ( j >i W+K 2 (d ( j }i K)W and is the inverse of the Kahler metric. This scalar potential 
holds for any two-derivative action. When (cfii) = for all i, the Weyl rescaling is unnecessary 
and we can determine the contributions to the scalar masses directly from the curvature coupling. 
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Let us illustrate this for the case of slow-roll inflation: we may compute the mass of a directly 
from the F-term potential (A. 3) to find m 2 = 2k 2 \J 7 x\ 2 = QH 2 . Alternatively, we can expand 
(A. 2) using 1Z = —12H 2 to find one contribution to the mass of size ^^.(1) = By expanding 

the exponential (A.l) to order (<3? + &) 2 X^ X, we find an additional contribution, m 2 ,^ = 2-ff 2 . 
Combining both terms, we find m 2 = + m^^) = 6-ff 2 , as before. 

Let us make a few comments on the role of the auxiliary fields M and b^: in general, they can 
be integrated out to produce Planck-suppressed couplings between X and $. For two-derivative 
actions, these effects have been included in writing the potential in eq. (A. 3). The contributions 
of these fields are only relevant when M or acquire vevs. However, as we explain in §2.4, these 
vevs are small (or even strictly zero) during inflation and give at most subleading corrections to 
the masses of fields. 

For two-derivative actions, the combined effects of curvature couplings, the auxiliary fields, 
and Planck-suppressed interactions, are all encoded in (A. 3). However, for higher-derivative 
theories, there is no compact formula that includes all these terms [25]. In the context of inflation, 
determining the scalar potential remains relatively simple, as the auxiliary fields are negligible. 
As a simple example, we studied theories with small sound speed, such as 

K Cs = l^ + &f\ Cl + --^d^d^} , (A.4) 

where ci = 1 + |(l — \\ and c 2 = j(l — -jjj). In [25] we derive the component Lagrangian for 
this theory 

M*\H\ 

where 

A = -\d^\ 2 + (in - ^\T X \ 2 ) o 2 , (A.6) 



C 2 = - (\d^\ 2 f - 2d ll ad u ad^8^ 

+ ^{( \n-^\Tx\ 2 ) [(d^c) 2 + (d^a) 2 } - \n„ v [d^ c d v ip c + d^ad v (x\ 

+ \ VjNv [d»y c d v y c + d»od v o] + d^U4> + \U\d^\ 2 ] . (A.7) 

Using <p c = (2M 2 l \H\) 1 ^ 2 t and = -2>H 2 g^ (de Sitter), we determine the mass of the partner 
of the inflaton, 

£c 3 = ■ ■ • + [ci [\n- H 2 ) +c 2 {^R- 2H 2 ) {d^tf + c 2 V^V u (dHd u t)] a 2 + ■ ■ ■ . (A.8) 
With V^VuidHd^) = 9H 2 , we find 

H 2 

m 2 = 6H 2 Cl - \2R 2 c 2 ~ x — + 6H 2 . (A.9) 

c i 

The cancellation of the leading term in (A.9) is completely accidental and not protected by any 
symmetry. The mass received contributions both from curvature and from the exponentiation of 
the rigid action. Neither the curvature couplings nor the form of the Kahler potential are fixed 
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in minimal supergravity. We could change either coefficient to eliminate the cancelation without 
breaking any symmetry of the action. For example, there are contributions to the mass from 
modifications to the Kahler potential via Planck-suppressed couplings of the form 



jd 2 QE(V 2 -8R)j^X^X K Cs d 3f3H 2 a 2 



ci + 



C2 



M 2 ,\H\ 



\d u 



(A.10) 



These contributions are model-dependent, but generically present. This contributes terms at the 
same order and change the numerical coefficient in eq. (A. 9) 



m 2 ~ 



H 2 



^>H 2 



(A.ll) 



Hence, generically small sound speed will lead to a parametrically enhanced mass for a. 

We refer the reader to our companion paper [25] for a more complete treatment of super- 
gravity of more general effective theories. 
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B Estimates for Quasi-Single-Field Inflation 



In this appendix, we will explain how to estimate both the amplitude of the bispectrum and 
the scaling of the squeezed limit for models of quasi-single-field inflation. In both cases, we will 
provide two different methods for estimating the results: 

i) We will use dimensional analysis to estimate the answer from the full bispectrum calculation. 
The advantage of this approach is that it essentially sets up the complete calculation and 
therefore there is little room for error in the estimate. 

ii) We will estimate the answer directly based on physical reasoning. The advantage of this 
approach is that it dramatically simplifies the process of estimating the results and provides 
a physical understanding of the origin of the effect. 

Of course, both approaches will yield the same results. 
B.l Amplitude of the Bispectrum 

We use the "in-in formalism" (see e.g. [52, 53]) to compute the bispectrum 



<C 3 (o)> = (n\ 



Texpfi f if/(r)dr) C/(°) Texp(-i f #/(r)dr) 

J -co J L J -co 



\n) , (B.i) 



where the subscript / indicates interaction picture fields and T denotes time ordering. We can 
calculate the bispectrum of any operator in perturbation theory by expanding the time-ordered 
exponentials to the appropriate order in the interaction Hamiltonian Hi. 

Direct estimation. As a concrete example, let us consider the original model of QSFI [11], 
with interaction J d^x^/^^g [ia 3 and mixing term J d^Xy/^^g m 3 Tra. The bispectrum generated by 
these interactions is schematically of the form 



(Cl) = Cfc(0) 3 [J / dr,a 3 m 3 7rV(7i) • J draW) , (B.2) 



44> 

i=l 

where ' denotes a derivative with respect to conformal time r. Here, we have ignored the details 
of the time ordering and the commutator structure that follows from eq. (B.l). We will be 
estimating the results by dimensional analysis of the integrals which will be insensitive to these 
subtleties — although these subtleties can be important to prove convergence of the integrals [11]. 
To estimate the integrals, we will need the mode functions for tt and a. Since the Goldstone 
boson is massless in the decoupling limit, its mode function in de Sitter space is 

vr fe (r) = -L= -(1 + ikr)e- ikT . (B.3) 

The mode functions for the massive field a are Hankel functions. When \kr\ <C 1, these simplify 
to 



.(_ T )3/2-v 



h^ =ff — where ^Vi-#- (B - 4) 



in 



2 
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This long-wavelength behavior will be sufficient for estimating the integrals, since the modes 
oscillate rapidly and suppress the integrals when \kr\ 3> 1. 

The first integral that we wish to estimate is the mixing term f dr a 3 m 3 7r / (T. Using eqs. (B.3) 
and (B.4), we can write this as 



-e - ikT 



I" dr Q 3 mV<x~ p P dr , = k ^ e~* k A ^ " ~ L f* dr k V*-» T -»-V\ 
J H J (-ffr) 3 V5P & Hj 

(B.5) 

As we stated above, when \kr\ S> 1, the integral is exponentially suppressed by the rapid os- 
cillations of the mode functions. Since the only scale in the problem is k, the integral will be 
dominated by r ~ A;" 1 . For large values of v one might worry that the integral is actually diver- 
gent as the upper limit of integration — > 0, but a more careful analysis shows that all power 
law divergences cancel [11]. In the following, we will take it as a given that there aren't any 
power law divergences. This ensures that the integral receives its dominant contribution from 
r ~ k^ 1 . Moreover, this observation also shows that these integrals are only weakly dependent 
on the upper limit of integration. The integral in (B.5) may then be estimated as 

dro 3 mV(T~^fc . (B.6) 
H 

This leads to the "Feynman rule" that each insertion of the mixing contributes a momentum- 
independent factor of p/H to the bispectrum amplitude, /nl- 

The interaction term, f dra 4 //<7 3 , can be estimated in a similar way, 

^W^f^t^ „ | /^-£^ . ,b„ 

We have again assumed that any power law divergences cancel and that the integral is dominated 
by t ~ /c" 1 . We see that each interaction vertex leads to a factor of p/H. 
Putting these estimates together, we find 

where = H 4 /(4M^\H\). This result is consistent with the complete calculation in [11]. 

Quick estimate. While the estimate that we just performed is robust, it relies on actually 
setting up the full calculation. We would like to have a faster and possibly more intuitive 
explanation for the result. In single-field models, we can often estimate the size of non-Gaussianity 

as 

/nl~ 7 /, (B.9) 

where the r.h.s is evaluated at horizon crossing. At first glance, it may not be clear how such an 
estimate would work for quasi-single-field inflation. First of all, the interaction £3 involves only 
a, while the power spectrum involves just ir. This problem is easily resolved by using the fact 
that the mixing term is independent of k. Specifically, the mixing term allows us to contract tt 
with a at the cost of a numerical suppression, p/H. In principle, the mixing term therefore allows 
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us to compute the ratio in (B.9). The second issue we face is that the a modes evolve outside 
the horizon. To make an estimate, we therefore need to specify the time at which the amplitude 
of a is measured. However, our previous estimates have shown that the dominant contributions 
to the integrals comes from the time of horizon crossing r ~ k^ 1 . Since /nl is defined as the 
amplitude of the bispectrum in the equilateral configuration (k\ = k2 = k%), we may evaluated 
each a at horizon crossing. This allow us to rewrite a in term of ir, 

°k~i/r -> ^(M^Hl) 1 / 2 7r k ^ 1/T = -|^vr^ 1/r . (B.IO) 

We can then estimate /nl from eq. (B.9), 

/NL &H A < ■ (B " n) 

We see that this estimate matches the more careful estimate in eq. (B.8). 

Generalizations. Above we reproduced the result for the specific QSFI example of [11]. In 
Section 4, we explored variations of QSFI, allowing for all possible interactions between n and 
a consistent with the symmetries of the effective theory of inflation. We now have the tools to 
understand the non-Gaussianity in all these examples and to recover the results of Table 1: 

First, we consider the alternative mixing term of the form 

a(M^\H\) 1/2 ird- . (B.12) 

By repeating the above analysis, we find that the integral associated with these mixing terms 
scales as ak°. The coupling a therefore now plays the role of p/H . The rest is unchanged. Hence, 
replacing p/H by a gives the result for the new mixing terms. 

Next, we consider models with different interactions, i.e. different couplings between tt and 
a and different derivative structures. Applying the methods of this appendix, we can determine 
Table 1. We find an enhancement in /nl by factors of A^ 1 for certain cubic interactions. This 
can be explained by the larger amplitude of a fluctuations at horizon crossing relative to the 
fluctuations of tt. From eq. (B.IO) we see that we gain by powers of A" 1 for every factor of a 
that appears in £3. For example, estimating /nl for the interaction fhad^ird^aa, we find 

, 1 £3 _ 1 mgd^TTd^aa fh a ir p 2 _ fh a p 2 m 1^ 

~ CC 2 ~ C M 2 yH (dj/ir) 2 C H 2 ~ H IP ' [ ' 6) 

As expected, /nl for this interaction is suppressed relative to (B.ll) by one power of A^. 
B.2 Squeezed Limit of the Bispectrum 

In the previous section, we showed how to estimate the amplitude of the bispectrum. We will 
now use a similar analysis to determine its scaling behavior in the squeezed limit. 

Direct estimation. The most reliable way to estimate the squeezed limit is by dimensional 
analysis. We again start with the schematic form of the bispectrum 

3 



(Ck 1 Ck 2 Ck 3 )(0) = <5(k 1 + k 2 + k 3 



CfciCfc 2 Cfc 3 (0) • \\ J ^i^^^'k^k, ■ j dTa A pa kl a k2 a ka 



(B.14) 
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where the delta function enforces k<i = k% in the squeezed limit, k± — > 0. As before, we will 
estimate the behavior of each piece separately: 

The easiest piece to understand is the contribution from the overall factor Cfci Cfc 2 Cfc 3 (0) , 
which scales as 

A 3 

CfciCfaCfcsCo) = ^4 • (B - 15) 

k 1 k 2 

Next, we consider the mixing terms. As we saw in the previous section, they scale as k°. 
In principle, they could therefore contribute to the squeezed limit through ratios like {kxjk^f 1 ■ 
However, this does not occur, since each mixing term is only a function of a single momentum. 
Specifically, the integral over the long-wavelength modes with momentum k\ takes the form 

draW ki a kl ~ | f dr^V'^/V** ~ | . (B.16) 

As before, we estimated the integral by taking the dominant contribution to the integral to come 
from r ~ k7 . The only way the scale &2 could have entered the estimate, is through an implicit 
^-dependence in the integration limit r*. However, as we explained above, the cancelation of 
divergences as r* — > 0, also implies that the integral is only weakly dependent on the limits of 
integration. As a result, any ^-dependence from the limit of integration is a subleading effect. 
The skeptical reader is encouraged to consult the full in-in calculation [11] to confirm this. 

The final contribution is from the interaction term. In the previous section, we found that 
this integral scales as k~ 3 / 2 , as required for scale invariance of the bispectrum. Repeating this 
estimate, but now being careful to distinguish the long mode, k%, from the short modes, k%, we 
find 



limj dr^a kl a k2 a k , ~ - J dr . (B.17) 

The integral contains two scales, k% and k2, so it may not be obvious which scale dominates. 
When |r| ~ k^ 1 S> k^ 1 , the short modes are still well inside the horizon and oscillate rapidly. 
Therefore, the integral is exponentially suppressed (this isn't obvious in eq. (B.17) because we 
have taken the mode functions in the limit where the modes are outside the horizon). It isn't 
until the short modes freeze out, I 7" | rs ~' A?2 j theft we get a significant contribution to the integral. 
Therefore, we can estimate the integral as if there were only one scale, &2, to get 

lim / dr ^a kl a k2 a k ~ ^ . (B.18) 

Putting all these estimates together, we find 

Km (OhCfaCfa) oc fc 3/2 + 49/2-, • (B - 19) 

This agrees with the result of [11]. 

Intuitive explanation. To gain some intuitive understanding for the scaling of the bispec- 
trum in the squeezed limit, it is useful to rewrite the result as 

So^w^U) ' (R20) 

local 
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The first term has the same squeezed limit as the local shape, /cj~ 3 . The second term modifies this 
scaling by the ratio (k\/k2) 3 ^ 2 ~ u ■ For u < 3/2 (or m a > 0), this leads to a suppression in the 
squeezed limit, while for v = 3/2 (or m a = 0) we recover the local shape. This suggests that the 
physics of the squeezed limit of QSFI is the same as the local shape, with a modification induced 
by the superhorizon behavior of the mode function for the massive field a. 

Let us try to understand this in a bit more detail. The squeezed limit of the local shape 
can be understood as a local modulation of the short-scale power due to the presence of the 
long-wavelength mode. For m 2 = in QSFI, we recover the same result. On the other hand, 
when m 2 > 0, the isocurvaton decays outside the horizon. Although the modulation of the 
power spectrum is still local in space, the amplitude of the long mode now depends on the time 
of horizon crossing. Specifically, the interaction plays an important role at horizon crossing of 
the short modes, k%. When the short modes cross the horizon, \k%T2\ ~ 1, the amplitude of the 
long mode, k\, is suppressed relative to its amplitude at its own horizon crossing, \k\T\\ ~ 1, 

M^) = Mn)(jM • (B.21) 

We hence understand that the deviation from the local shape in (B.20) is the result of the decay 
of the long-wavelength mode between the time it crosses the horizon and the time the short 
modes cross the horizon. The mixing term also plays an important role in this interpretation. 
It allows the conversion of the massive mode a into the massless mode £ = —Htt shortly after 
horizon crossing. If one simply computed the bispectrum for a using the cubic interaction, one 
wouldn't find the same squeezed limit. Unless the a fluctuations are converted into a massless 
mode at horizon crossing, they will continue to decay outside the horizon. This would lead to 
additional factors of k\ and k<i related to the decay of the amplitude from horizon crossing until 
the end of inflation. 

Generalizations. The behavior in the squeezed limit is robust to changes of the model: 
The specific form of the mixing term is unimportant, as its only role in life is to convert 
a massive mode into a massless mode. As we discussed in the previous section, even when we 
include more derivatives in the mixing term (e.g. 7r<r), the integral still scales like k°. As a result, 
the estimate of the mixing integral is the same in all cases and does affect the squeezed limit. 
Changes to the interaction also do not alter the squeezed limit, as long as there is at least one 
field with no derivatives acting on it. The behavior in the squeezed limit, is the result of a "local" 
correlation of the long and short modes. When there are no derivatives acting on the field, the 
short modes are sensitive only to the amplitude of the long mode, at the time of horizon crossing. 
In the case, where all fields are acted on by derivatives, the short modes are only sensitive to 
gradients of the long mode, and therefore the squeezed limit is further suppressed by a factor 
of k\. 
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